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SUMMARY 

Axial segmentation of acoustically absorbing liners in rectangular, 
circular or annular duct configurations has proven to be a very useful 
concept to obtain higher noise attenuation with respect to the bandwidth of 
absorption as well as the maximum attenuation. As a consequence, advanced 
liner concepts have been proposed which induce a modal energy transfer 
in both cross-sectional directions to further reduce the noise radiated 
from turbofan engines. However, these advanced liner concepts require 
three-dimensional geometries which are difficult to treat theoretically. 

In the present paper, a very simple three-dimensional problem is investi- 
gated analytically. The results show a strong dependence on the positioning 
of the liner for some incident source modes while the effect of three- 
dimensional segmentation appears to be negligible over the frequency range 
considered. 


t This research was conducted while the author held a National Research 
Council Senior Research Associateship at the NASA Langley Research 
Center, Hampton, Virginia. 
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INTRODUCTION 

In order to meet current and future noise regulations of commercial 
jet transport aircraft, considerable research effort has been expended over 
the past years in an attempt to attenuate the noise radiated from turbofan 
engines. With most design possibilities for minimizing the noise generated 
at the source being exhausted, acoustic lining became one of the more 
established approaches to further reduce the noise. Advanced duct liner 
concepts have been developed in order to achieve more absorption with less 
weight and volume of acoustic treatment starting with the concept of axially 
segmented liners originally introduced by Lansing and Zorumski [1]. 

Several theoretical and experimental investigations have proved that multi- 
segmented liners can be superior to uniform liners with regard to the band- 
width of absorption as well as the maximum attenuation. The basic idea 
behind the concept is that the impedance discontinuities of the various 
segments are used to redistribute the acoustic energy into different, 
usually higher order} acoustic modes which are absorbed more efficiently. Of 
course, the results are strongly dependent upon the given sound source. 

In [1], as well as most other papers treating multisegmented liners 
(see for example [2], [3], [4]), either two-dimensional or annular duct 
geometries have been employed which resulted in one mode number remaining 
constant. From the fact that, in an annular lined duct, the attenuation 
is higher not only for higher order radial but also for higher order 
circumferential modes, it is only natural to expect higher attenuation if 
the modal structure is broken up in both mode numbers since then several 
more modes are available for an energy transfer. The use of circumferentially 
as well as axially segmented liners leading ultimately to the concept of 
liners with continuously varying impedance in all directions is based upon 



this idea. Aside from manufacturing difficulties, the main problem with 
these advanced liner concepts is that they involve three-dimensional 
geometries which are difficult to treat theoretically. 
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The present investigation is a first attempt to examine the effects 
of three-dimensional segmentation analytically. By looking at a very simple 
three-dimensional problem, namely two axially segmented uniform liner pairs 
which are tilted by 90 degrees with respect to each other in an otherwise 
infinitely long, hard walled rectangular duct, the effect of modal break up 
in both cross-sectional directions can be studied. At the same time, this 
problem constitutes an excellent testing ground for the analytical approach 
of combining the mode matching and Wiener-Hopf technique. This method is 
very useful in solving certain three-dimensional problems by matching essentially 
two-dimensional solutions. Disappointingly, the few sample results that were 
evaluated numerically show practically no increase in sound attenuation due 
to three-dimensional segmentation. Two-dimensional segmentation appears to 
be equal or even superior in the examples considered as long as the liners 
are positioned correctly. Thus the results reemphasize the fact that, for 
a given incident mode, the sound attenuation can be increased substantially 
if the liner is positioned such that the wave front normal points into the liner. 

GOVERNING EQUATIONS 

The geometry of the problem is schematically depicted in Figure 1. In 
an infinitely long, hard walled rectangular duct of constant cross-section, 

WxH, two liner pairs, separated by a distance AL = - L-|, are mounted on 

the side walls and are tilted by 90 degrees with respect to each other. A 
uniform axial mean flow of Mach number M = W^/a^ >0 is assumed in the positive 
x - direction. Viscosity, thermal conductivity and all nonlinearities are 



ignored. Introducing a dimensionless perturbation velocity potential $ 0 
such that 
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the governing convected wave equation is 
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with the corresponding expression for the perturbation pressure p 
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L re ^ is a reference length (in our case, the duct height H) while the reference 
velocity U ^ is taken to be the uniform mean flow velocity W^, if mean 
flow is present, or the ambient sound speed a^ if the undisturbed medium 
is at rest. The sound source is assumed to have harmonic time dependence 
exp (iwt). 

After applying a Prandtl-Glauert type transformation and simultaneously 
nondimensionalizing the spatial coordinates 
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it is advantageous to introduce the amplitude function cj>(x,y,z) defined by 
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This reduces the governing equation (2) to the Helmholtz equation 
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which has to be solved for the axially segmented liner configurations shown 
in Figure 2 in nondimensional ized notation, i.e. h = H/L rgf , w = W/L ref 
and i = L/(L ref V 1 - M 2 ) . In Figure 2a, the two liner pairs are mounted 
on the same side walls while the remaining two walls are hard. This is a 
two-dimensional configuration since the mode number in the y-direction does 
not change throughout the whole lined section. In Figure 2b, the second liner 
pair is tilted by 90 degrees. It is this three-dimensional configuration 
we are mainly interested in here. 

In the hard walled sections I, II and III (see Figure 2) the sound field 
can be described by the sum of acoustic modes n (y, z). The method of 
separation of variables together with the application of the hard wall boundary 
condition of zero normal velocity leads to (cf. [2]) 


x 4 0: 


(I) 


(m,i)=E (V l 5 a ( '> e' lY ™." x V (y, 2), 

r,s r,s m=Q n=Q m,n m,n 


a, £ x $ H 0 : 

I v x £ 00 co 

* (II) (x.y. z| - 5 J e^. n T (y, z) + 

m=0 n=0 


CO 00 


* } ) A (2) e -lYm,n ( l 2 ‘ x > T (y, z) 

io S=0 m ’ n 


x £ £ 2 : 


00 00 


(,n, (*.y. Z )-}„ L e'*.n <* - h> Vn (y.*) + 


m=0 n=0 


+ E e 
r,s 


■■*Yr,s ( x " ^3) 



G 


Here the hard wall duct eigenfunctions and eigenvalues are 
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Assuming only a single incident duct wave (due to the linearity, the 

solutions can always be superimposed), E , E denote the modal amplitudes 

r,s r*s 

of the prescribed incoming hard wall duct waves impinging from upstream and 

downstream respectively. A ^ , A ^ , and B ^ , B ^ are the modal 

m,n m,n m,n m,n 

amplitudes of the waves in front of and behind the v-th lined segment, 
propagating upstream or downstream (see Figure 2). 
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are the corresponding modal propagation constants. For K <K^ ^ , the duct 

wave is attenuated (cut off) while for K > K the duct wave propagates. In 

rn 5 ri 

the latter case one has to distinguish two cases (cf. [2]): 

a) K > n / V~T-M 2 : The corresponding phase velocities have different 

signs for the upstream and downstream propagating wave. 

b) K < K <K / VI- M 2 : The two phase velocities have the same 

' m,n m,n * 

sign and energy considerations are necessary to determine the direction of wave 
propagation (see [5] or [6], p. 161). 

For the lined sections, a locally reacting lining material of specific normal 
acoustic impedance ^ is assumed. For our numerical examples, the crossover 
frequency impedance model 



7 


^ = R* (1 + i f*/ f 0 J " 1 ctg ( 2n f * d *^ < 5 ) 

for honeycomb liners with an impervious backing and a thin sheet of absorptive 

facing material is employed. R* and f 0 * = f 0 L re f/a ro are the specific 

acoustic resistance and the nondimensional characteristic frequency (crossover 

frequency) of the facing sheet while d* = d/L ref is the dimensionless depth 

of the honeycomb structure, f* = f L ^./a^ denotes the frequency parameter. 

Since we allow a uniform mean flow, the medium inside and outside the 

liner is assumed to be separated by an infinitely thin vortex sheet of finite 

axial length across which the pressure remains continuous. Furthermore, the 

kinematic condition of equal particle displacement on both sides of the vortex 

sheet must be imposed. The solution of the problem of sound transmission through 

a finite length acoustically treated duct section is assumed to be known in terms 

of the various reflection and transmission factors R and T 

m,n;y,v m,n;u,v. 

The details of the solution are unimportant as far as the present work is con- 
cerned to the extent that the finite length liner can be taken as a "black box" 
with known reflection and transmission characteristics. 

GENERAL FORMULATION IN TERMS OF THE MODAL REFLECTION AND TRANSMISSION FACTORS 
The concept of reflection and transmission factors has proven very useful 
in several engineering disciplines. In acoustics, most muffler design procedures 
are based on the one-dimensional transmission-line approximation and a very 
nice introduction to multidimensional problems in ducts was given by Zorumski [7 J . 
The reflection and transmission factors relate the modal amplitudes in front 
of and behind a nonumform duct section. This nonuniform duct section may be a 
discontinuity in the physical properties, like a change in duct cross section or 
wall admittance, or may be of finite axial extent. In this paper, we are mainly 
concerned with the latter. Referring to a simple liner pair in Figure 2, 
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R , for example, denotes the complex amplitude of the (m,n) th 

m,n;y,v ^ r K 

reflected wave due to a unit amplitude (y,v) th wave impinging from upstream. 

Similarly T denotes the complex amplitude of the (m,n) th transmitted 

hi > n *y 3 v 

wave due to a unit amplitude (y,v)th wave impinging from upstream. Similar 
definitions apply for the incident wave coming from downstream indicated by a 
left pointing arrow in the superscript. 

If all four walls of the two lined sections in Figure 2 are treated with 
acoustic material, one can write for the first lined section 
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and similarly for the second liner segment 
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With zero mean flow the reflection (and transmission) factors R . and 

m 9 n 9 y 9 v 

R ' v ' are equal. Throughout the intermediate hard walled section II, each 
m,n;y,v ^ 3 

mode propagates independently such that 


A (2) g (1) e Tym,n U 2 - ) 

m 5 n m 5 n 


( 8 a) 


B ( 1 ) = a ( 2 ) e " iYm ’ n (£ 2 - *l) • ( 8 b) 

m,n m,n 

Equations ( 8 ) can be used to eliminate A anc j $ (P in ( 6 ) and (7). Then 

m,n m,n 

( 6 b) and (7a) constitute a two-fold infinite set of determination equations for 
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the two-fold set of infinitely many unknown amplitudes B ' ' and A ' 

m,n m,n 



9 


m,n = 0, 1 , 2, . . . , 
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This infinite set of equations can be truncated a few modes beyond the modal 

cut off numbers. Therefore, once the reflection and transmission factors are 

known and the source characteristics are prescribed in terms of the incoming 

modal amplitudes E and E , the solution of the truncated set of 

r 5 s r,s 

equations (9) is straight forward. 

The computation of the reflection and transmission factors constitutes 
the main problem. Unfortuately , to my knowledge there exist no analytical results 
for the reflection and transmission factors of three-dimensional rectangular 
ducts lined uniformly over a finite axial length on all four walls although the 
mode matching method should provide one means to do so. With regard to numerical 
techniques, for example, Watson's method [8] could be employed. However, in 
several applications only two opposing side walls are lined such that one mode 
number can be held fixed while computing the corresponding reflection and 
transmission factors. This is the case for the examples depicted in Figure 2. 

Thus the solution procedure will now be outlined for the special case of two tilted 
liner pairs. 

TWO LINER PAIRS TILTED BY 90 DEGREES WITH RESPECT TO EACH OTHER 
In the first lined section of Figure 2b, the walls are hard at y = 0 and 

-*■ +■ (i) 

y = w. Hence, for a given incoming mode E or B , the mode number r or 

r )S in 9 n 
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m in the y-di recti on will be unaffected by the presence of the first liner and 

will equal that of the incoming mode. For the quasi-two-dimensional problem 

shown in Figure 2a, this is true all the way through both lined sections such 

that r is fixed by the incoming mode. For the second lined section of Figure 2b, 

-*■ (?) 

the walls at z = 0 and z = h are hard and for a given incoming mode, A v ' 

J J m,n 

4" 

or E , the mode number s or n in z-di recti on will equal that of the 
i >s 

incoming mode. The three-dimensional effect is therefore due to the interaction 

of the two liner pairs. In the intermediate hard walled section II, an incoming 

"*■ (?) (2) 

mode A produces reflected amplitudes n with n = const., while an 

( i ) * ->■ ( i ) 

incoming mode B v ' produces reflected amplitudes B ' ' with m = const. 

3 m,v K m,n 

All reflection and transmission factors of the first lined section can be 
computed by keeping the mode number m constant and equal to that of the 
corresponding incident mode; and all reflection and transmission factors of the 
second liner section can be computed by keeping the mode number n constant 
and equal to that of the relevant incident wave. 

Equation (6) can therefore be rewritten in the form 
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Here 6 denotes Kronecker's delta function. 
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Again employing the uniform duct transmission equatiors (3) , the infinite 
system (9) now reduces to the special form 
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r,s are the prescribed mode numbers of the incident waves. Using matrix 
partitioning this set of equations may be written 
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NUMERICAL RESULTS 

The quantities of most practical interest are the reflection and transmission 
coefficients . These are the ratio of the reflected or transmitted power in the 
(m,n)th cut on mode over the incident power in the (r,s)th cut on mode. For 
an (r,s) duct wave incident from upstream, they are defined by 
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m,n;r,s = 
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Similar definitions apply for a mode incident from downstream. The total 
reflected or transmitted power for a single incoming mode (r,s) is then 
obtained by simply summing the contributions of all cut on (m,n) modes, i.e. 


r,s 


- I 


m,n/ 


m,n;r,s 


cut on 


r,s 


ran/ 


m,n;r,s 
cut on 


The acoustic effectivity of a liner combination is usually measured by the 
power attenuation or insertion loss AP. For a single incoming mode (r,s), 
it is defined by 

AP [dB] = 10 log (1/T r>s ) . 

In our numerical examples, the mean flow Mach number is taken to be zero 
such that the relevant reflection and transmission factors for a simple liner 
pair can be computed by the method outlined in [4]. Originally it was intended 
to compute the reflection and transmission factors for uniform mean flow 
("plug" flow) by extending the Wiener-Hopf method described in L4]. It turned 
out that the extension is not trivial leading to the appearance of instability 
waves hitherto neglected in all mode matching approaches. Therefore a more 
thorough study of the mean flow problem is necessary which will be reported 
separately. 
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The infinite set of equations (12) is evaluated numerically by truncating 
it, in general, three modes beyond the corresponding cut off mode. Valuable 
computational checks are provided b w the requirements that the total transmitted, 
reflected and absorbed powers have to equal the incident power, i.e. 


as well as by the reciprocity theorem in the form 
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This reciprocity theorem is valid only for M = 0 but applies to cut off 
modes as well as to cut on modes. The relevant absorption coefficient a 

r ,s 

can be computed from the solution in the lined section (see [4]). 

The numerical examples are chosen such as to provide a comparison with 

previous results [4] whenever possible. In particular, the two liners in 

each pair are assumed to be identical and to have the same axial length 
2 1 7 

Ll/H = — • To avoid geometric effects, the duct is assumed to be of square 

cross section H x H, H = 18.73 cm. The crossover frequency impedance model (5) 
with f * = 8.14 is used for all examples. 

In Figures 3 to 5 the power attenuation is plotted as a function of frequency. 
First, to check the procedure, the two-dimensional configuration of Figure 2a is 
evaluated. Assuming an incident fundamental mode (0,0), the result is shown as 
curve A in Figure 3 and is identical to the one given in [4], From a computational 
point of view, the present approach, using mode matching in conjunction with the 
3-part Wiener-Hopf technique, is preferential to the n-part Wiener-Hopf approach 
of [4] because one and the same Wiener-Hopf subroutine can be used for an arbitrary 
number of lined sections and only the fairly straight forward matrix equation (12) 
has to be reformulated if more liner sections are added. Also shown in Figure 3 
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are the results for the case when the two liner pairs are separated by a hard 
walled section of length L] (curve B). Only minor differences are observed 
due to the interference of waves in the intermediate section. Finally, the 
second liner pair is tilted by 90 degrees (curve C) and again only minor changes 
occur. As a matter of fact the power attenuation is even slightly below the 
one for two-dimensional segmentation over most of the frequency range shown. 

In Figure 4, the incident modes (1,0) and (0,1) are considered. Curves A 
and B show the results for two-dimensional segmentation. Depending on whether 
the normal wave front vector points into the lined wall (curve B) or is merely 
tangent (curve A), substantial differences are observed in the power attenuation 
with a corresponding shift in the tuning frequencies of the two liner segments. 

Curves C and D show the solution for three-dimensional segmentation, i.e. the 
second liner pair is tilted. It appears that the effect of correct positioning 
of the liners by far outweighs the effect of three-dimensional segmentation. Up 
to the second resonance of the second liner segment at f* ~ 2, the curves for 
the tilted configuration essentically follow the ones for the corresponding sinqle 
two-dimensional liner section and then remain at the higher level. Translated 
to the situation of an annular turbomachinery duct, where mainly circumferential 
modes are excited, this means that radial liners would be vastly more efficient 
acoustically than annular liners of the same absorbing area. In a different context, 
the advantages of the concept of radial liners were pointed out by Abdelhamid 
[9]. 

Finally, in Figure 5, the incoming mode is taken to be (1,1) and hence a 
different liner combination with tuning frequencies above f* = 0.7 had to be 
chosen. A summary of the cut on frequencies of the various modes for all 
cases is presented in Table- 1. For comparison, the sound attenuation of the 
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fundamental mode is also included (curve A). Curves B and D give the results 
for two- and three-dimensional segmentation respectively, for the case with the 
better impedance match, i.e. the thinner liner is closer to the source. Curves 
C and E depict the results for the corresponding cases with the two liner 
segments interchanged. It is somewhat surprising that, even here, the two- 
dimensionally segmented duct outperforms the duct with three-dimensional 
segmentation over most of the frequency range, similar to the situation in Figure 3. 

CONCLUDING REMARKS 

The simple three-dimensional problem of two axially segmented liner 
pairs tilted by 90 degrees with respect to each other in a rectangular duct 
was used to study the effect of three-dimensional acoustic liner segmentation. 

Only a few sample results have been evaluated. However, none of these indicate 
that the effect of modal energy transfer in both cross-sectional directions can 
be used to obtain a substantially higher sound attenuation per given liner 
length. The effect of correct liner positioning, i.e. such that the wave normal 
points into the liner and is not merely tangential, appears to outweigh the 
effect of three-dimensional segmentation by far. For the situation of a turbofan 
exhaust duct, where the sound generation by blade rows leads to the excitation 
of mainly circumferential modes, this would mean that radial liners are 
acoustically superior to annular liners. 

Although the results with regard to three-dimensional segmentation are very 
disappointing, the analytical method employed has more far reaching 
applications already envisaged by Zorumski [7]. In the present paper, the 
mode matching method is used only to link different finite length sections 
via the respective transmission and reflection factors. Very often the 
geometry of those finite length sections is such that only an interaction with one 
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mode number can occur. In such a case, the various reflection and transmission 
factors can be computed by two-dimensional methods like the Wiener-Hopf technique. 
Linking these by means of the mode matching method makes it possible to build 
up three-dimensional problems by means of two-dimensional solutions. One such 
example is the interaction of sound generated in a blade row and sound reflected 
and absorbed by an adjacent liner. Other applications include the interaction 
of different blade rows (or stages), not only with respect to their acoustics 
like rotor shielding in turbofan engines or sound transmission through several 
compressor stages, but also with respect to their aeroelastic behavior since 
blade vibrations are also communicated via acoustic modes. 
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Latin symbols 

A (V) 


m,n 


m,n 

d* = d/L 


ref 


'r,s 


f*=fL ,/a 
f o* -f o L ref^ a «> 


H 

h - H/L ref 

i-VT 

i 


K = “ L ref /(a ~ V 1 - M 2 ) 
!<,„ = m n/w 


K n = n n/h 


4 


,n V^m 


2 2 
+ K n 


L r 4’ 4 


ref 

1 = h/O-ref V 1 “ 

M = W /a_ 


P 

P ro 

R* 

R 

m,n;y,v 


LIST OF SYMBOLS 

modal amplitude upstream of v-th liner segment 
ambient speed of sound 

modal amplitude downstream of v-th liner segment 
dimensionless honeycomb depth 
modal amplitude of incident wave 
frequency 

frequency parameter 

reduced characteristic frequency of liner facing 
sheet 

duct height 

dimensionless duct height 

imaginary unit 

unit vector in x-di recti on 

reduced circular frequency parameter 

m-th hard wall eigenvalue in y-direction, m = 0,1,2,... 

n-th hard wall eigenvalue in z-direction, n = 0,1,2,... 

axial lengths; see Figure 1 
reference length (= H) 
reduced axial length L-| ; see Figure 2 
mean flow Mach number 
dimensionless acoustic pressure 
ambient pressure 

specific acoustic resistance of liner facing sheet 

(m,n)th modal reflection factor due to (y,v)th incident 
mode 
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Re ( z ) 
T m,n;p,v 


ref 


W 

U 


w = W/L ref 
X, Y, Z 


x, y, z 


real part of z 

(m,n)th modal transmission factor due to (y,v)th 
incident mode 

time 

reference velocity; U ref = a ~ if M U ref = 
W if MM 

CO 

dimensionless acoustic velocity vector; see 
equation (1) 

duct width 

uniform mean flow velocity 
dimensionless duct width 
rectangular coordinates 

reduced spatial coordinates; see equation (3) 


Greek symbols 


/s 



V K 2 - K 2 , K > K 
’ m,n m,n 

- K? • K < K m,n 


(m,n)th eigenvalue of hard 
walled rectangular duct 


AP = 10 log (l/T rjS ) 

6 

m,n 





V 


p 


r,s 


p m,n;y,v 


x 


r,s 


T 


m,n;y,v 


power transmission loss for (r,s)th mode incident 

Kronecker's delta; has a value of one if m = n and 
of zero if m / n 

specific acoustic normal wall impedance; see equation (5) 
specific heat ratio 

total reflection coefficient for (r,s)th incident 
mode 

modal reflection coefficent for (y,v)th incident 
mode 

total transmission coefficient for (r,s)th incident 
mode 


modal transmission coefficient for (y,v)th incident 
mode 
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$ (X, Y, Z, t) dimensional velocity potential; see equation 

cp(x, y, z) amplitude function; see equation ( 4 ) 

„(y> z ) = cos (Krj) cos (K_z) (m,n)th eigenfunction of hard walled rectangular 
in 5 n nr n . , 

duct 

to = 2nf circular frequency 



quantity referring to wave propagating in the 
downstream direction 


C) 


quantity referring to wave propagating in the 
upstream direction 
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LEGENDS FOR ILLUSTRATIONS 


Figure 1.- Schematic view of two axially segmented duct liner pairs tilted 
by 90 degrees with respect to each other. 


Figure 2.- Geometry of axially segmented liner in nondimensional notion: 

(a) two-dimensional segmentation; (b) three-dimensional segmentation 

Figure 3.- Power attenuation versus frequency with (0,0) mode incident. 

Figure 4.- Power attenuation versus frequency with (1,0) or (0,1) mode incident. 
(Liner data as in Figure 3). 


Figure 5.- Power attenuation versus frequency with (0,0) mode, (Curve A), or 
(1,1) mode incident. 


Legend for Table 


Table 1.- Cut on frequencies f* for various modes (m,n) in a square duct 
(zero mean flow) . 



m 

0 

1 

2 

3 

4 

1 

5 

0 

0 



1.5 

2.0 

2.5 

1 

0.5 

0.7071 


1.5811 

2.0616 

2.5495 

2 

1.0 

1.1180 

1.414 

1.8028 

2.2361 

2.6926 

3 

1.5 

1.5811 

1.8028 

2.1213 

2.5 

2.9155 

4 

2.0 

2.0616 

2.2361 

2.5 

2.8284 

3.2016 

5 

2.5 

2.5495 

2.6926 

2.9155 

3.2016 

3.5355 


Table 1.- Cut on frequencies f* for various modes (m,n) in a 
square duct (zero mean flow). 
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